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on Entropic Dynamics (ED) is proposed. It is shown that the hyperbolicity of a non-maximally 
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an arbitrary system of 3N degrees of freedom leads to linear information-geometric entropy growth 
and to exponential divergence of the Jacobi vector field intensity, quantum and classical features of 
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The unification of classical theory of gravity with quantum theories of electromagnetic, weak and strong forces is 
f""^ , one of the major problems in modern physics. Entropic Dynamics (ED) namely the combination of principles 
of inductive inference (Maximum relative Entropy Methods, ME methods) 0, 0, 0] and methods of Information 
^ I Geometry (Riemannian geometry applied to probability theory, IG) Q, is a theoretical framework constructed on 
. 1 statistical manifolds and it is developed to investigate the possibility that laws of physics, either classical or quantum, 
might reflect laws of inference rather than laws of nature. Examples of dynamics that can be deduced from principles 
of probable inference are not absent in physics. The theory of thermodynamics [6[ and to a certain degree, quantum 
mechanics [7[, are examples of fundamental physical theories that could be derived from general principles of inference. 
In constructing an ED-model, the first step is to identify the appropriate variables (relevant information) describing the 
system and thus the corresponding space of macrostates. This is the most delicate step because there is no systematic 
', way to search for the right macro variables; it is a matter of intuition, trial and error. Once the information and 
normalization constraints are identified, using ME methods, the probability distribution characterizing the system 
£N) ' can be computed. Finally, using IG methods, a Fisher-Rao information metric [1, [9[ can be assigned to the space 
of macrostates of the system. Given the Fisher-Rao information metric, the geometric structure of the manifold 
underlying the ED can be studied in detail: metric tensor, Christoffel connection coefficients, Ricci and Riemann 
statistical curvature tensors, sectional and Ricci scalar curvatures, Jacobi and Killing fields can be calculated. ME 
methods are inductive inference tools. They are used for updating from a prior to a posterior distribution when new 
information in the form of constraints becomes available. Basically, information is processed using ME methods in the 
framework of Information Geometry. The ED model follows from an assumption about what information is relevant 
to predict the evolution of the system. In this work, we focus only on reversible aspects of the ED model. In this case, 
given a known initial macrostate and that the system evolves to a final known macrostate, we investigate the possible 
trajectories of the system. Given two probability distributions, a notion of "distance" between them is provided by 
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In this paper, our objective is to report some relevant results obtained in the realm of chaos theory (hypersensitivity 
to initial conditions) using the ED formalism. It is known there is not a well defined unifying characterization of 
chaos in classical and quantum physics [Hi [HI]. In the Riemannian Geometric Approach [T3, 1X3] to classical chaos, 
the search for a link between the Jacobi field intensity and the Ricci (sectional) curvature of the dynamical manifold 
is under investigation [Til ]. In the Zurek-Paz criterion of quantum chaos [l5| . instead, the search for a potential link 
between the linearity of the entropy growth and the curvature of the dynamical manifold underlying chaotic systems is 
still open. In our information-geometric approach, it is shown that these three indicators of chaos (Curvature- Jacobi 
Field Intensity-Entropy) are linked: the hyperbolicity of a 6iV-dimensional statistical manifold Ai s underlying an ED 
Gaussian model leads to linear information-geometric entropy growth and to exponential divergence of the Jacobi 
vector field intensity, quantum and classical features of chaos respectively. We assume that the arbitrary (physical- 
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biological) system under investigation has 3N degrees of freedom \ Xa \ , each one described by two pieces 

L J o=l, 2, 3 

of relevant information, its expectation value (x^\ and variance Ax^ = f J I (x^ — (x^^ Y This leads to 




consider an ED model on a 6-/V-dimensional statistical manifold M. s . First, we show that M. s has a constant negative 
Ricci curvature proportional to the number of degrees of freedom of the system, 1Zm 3 = — 3iV. Second, we suggest 
the information-geometric analog of the Zurek-Paz quantum chaos criterion. It is shown that the system explores 
statistical volume elements onM s at an exponential rate. We define an information-geometric entropy (IGE) of the 
system, Sm 3 - We show that <Sm s increases linearly in time (statistical evolution parameter) and it is proportional to 
the number of degrees of freedom of the system and to the information-geometric analogue of the Lyapunov exponents 
Finally, we show that the geodesies on the manifold A4 S arc described by hyperbolic trajectories. Using the 
Jacobi-Levi-Civita (JLC) equation for geodesic spread, it is shown that the intensity of the Jacobi vector field intensity 
J M s diverges exponentially (standard feature of classical chaos) and it is proportional to the number of degrees of 
freedom of the system. In conclusion, the Ricci scalar curvature 7Zm s , the information-geometric entropy Sm, and the 
Jacobi vector field intensity Jm s are proportional to the number of Gaussian-distributed microstates of the system. 
The relevance of this proportionality will be discussed in some detail in section III of this article. 

The layout of this paper is as follows. In section II, we describe the ED Gaussian model being studied. In section 
III, we introduce the information-geometric indicators of chaos for our theoretical model. Finally, in section IV we 
present our final remarks. 



II. ENTROPIC DYNAMICAL GAUSSIAN MODEL 



We consider an ED model whose microstates span a 3iV-dimensional space labelled by the variables 
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with a = 1, 2, 3 . We assume the 
and the variance 



only testable information pertaining to the quantities x^ consists of the expectation values I x^ 

Ax^K The set of these expected values define the 6JV-dimensional space of macrostates of the system. A measure 
of distinguishability among the macrostates of the ED model is achieved by assigning a probability distribution 

P(x e) to each 67V-dimensional macrostate 6 d = j { a a) ^ 6^ j = {((^ Q) ), A ^ Q) )} 3Ar with 

a = 1, 2,...., N and a = 1, 2, 3. The process of assigning a probability distribution to each state provides Ms with 
a metric structure. Specifically, the Fisher-Rao information metric defined in is a measure of distinguishability 
among macrostates. It assigns an IG to the space of states. Consider an arbitrary physical system evolving over a 

37V-dimensional space. The variables } label the 37V-dimensional space of microstates of 

the system. We assume that all information relevant to the dynamical evolution of the system is contained in the 
probability distributions. For this reason, no other information is required. Each macrostate may be thought as a 
point of a 6iV-dimensional statistical manifold with coordinates given by the numerical values of the expectations 
(^O^ and ( 2 )#^. The available information can be written in the form of the following 6N information constraint 
equations, 
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The probability distributions Pa in |T]) are constrained by the conditions of normalization, 
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Information theory identifies the Gaussian distribution as the maximum entropy distribution if only the expectation 
value and the variance are known [l7| • ME methods @, d, Q allow us to associate a probability distribution P (jl ^ 
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to each point in the space of states O. The distribution that best reflects the information contained in the prior 
distribution m (x) updated by the information (\Xa^\ , Axi™^ is obtained by maximizing the relative entropy 
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where m[X) is the prior probability distribution. As a working hypothesis, the prior m fxj is set to be uniform since 

we assume the lack of prior available information about the system (postulate of equal a priori probabilities). Upon 
maximizing (|3]), given the constraints (fTJ) and ([2|), we obtain 
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and, in standard notation for Gaussians. 
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. The probability 

distribution (|U) encodes the available information concerning the system. Note that we have assumed uncoupled 

constraints among microvariables x£*\ In other words, we assumed that information about correlations between the 
microvariables need not to be tracked. This assumption leads to the simplified product rule (|4]). However, coupled 
constraints would lead to a generalized product rule in Q and to a metric tensor ([7]) with non-trivial off-diagonal 
elements (covariance terms). Correlation terms may be fictitious. They may arise for instance from coordinate 
transformations. On the other hand, correlations may arise from external fields in which the system is immersed. In 
such situations, correlations among effectively describe interaction between the microvariables and the external 
fields. Such generalizations would require more delicate analysis. 

We cannot determine the evolution of microstates of the system since the available information is insufficient. Not 
only is the information available insufficient but we also do not know the equation of motion. In fact there is no 

(a) 



standard " equation of motion" . Instead we can ask: how close are the two total distributions with parameters (/j} a 
a^) and (pa + dpi" , 0a°^ + d<7a )? Once the states of the system have been defined, the next step concerns 
the problem of quantifying the notion of change from the macrostate O to the macrostate O + d<d. A convenient 
measure of change is distance. The measure we seek is given by the dimensionless "distance" ds between P (^X O^ 

and P(X e + d(=>), 



ds 2 = g^dBW 



where 



dXP 



8\ogP(x e)dlogP(x e) 



(6) 
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36^ 8Q» 

is the Fisher-Rao metric [1,0. Substituting ([4]) into J7J, the metric onX s becomes a QN x 67V matrix M made 
up of 3iV blocks M2x2 with dimension 2x2 given by, 
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with a = 1, 2,...., N and a = 1, 2, 3. From Q, the "length" element © reads, 
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We bring attention to the fact that the metric structure of Ai s is an emergent (not fundamental) structure. It arises 



only after assigning a probability distribution P I X 



] to each state 0. 



III. INFORMATION-GEOMETRIC INDICATORS OF CHAOS 

In this section, we introduce the relevant indicators of chaoticity within our theoretical formalism. They are the 
Ricci scalar curvature 71m s ( or i more correctly, the sectional curvature K.m s H)j the Jacobi vector field intensity 
J Ms an d the IGE Sm s - 



A. Ricci Scalar Curvature 



Given the Fisher-Rao information metric, we use standard differential geometry methods applied to the space of 
probability distributions to characterize the geometric properties of M. s . Recall that the Ricci scalar curvature 1Z is 
given by 

n = g» v R pv , (10) 
where g^ v g vp — S p so that = (g^)^ 1 - The Ricci tensor is given by, 

R/iv = 9 e T^ v — d u T ^ £ + r^ v T eTI — r^r^^. (11) 

The Christoffel symbols appearing in the Ricci tensor are defined in the standard way, 

r£„ = \g pe {d^g ev + d v g pB - d e9pu ) . (12) 
Using §§§ and the definitions given above, we can show that the Ricci scalar curvature becomes 

Km s = -3JV < 0. (13) 

From (|13[) we conclude that M. s is a 6iV-dimensional statistical manifold of constant negative Ricci scalar curvature. 
A detailed analysis on the calculation of Christoffel connection coefficients using the ED formalism can be found in 
[3 H3|- Furthermore, it can be shown that Ai s is not a pseudosphere (maximally symmetric manifold) since its 
sectional curvature is not constant. As a final remark, we emphasize that the negativity of the Ricci scalar TZms 
implies the existence of expanding directions in the configuration space manifold M s . Indeed, since 1Zm s i s the sum 
of all sectional curvatures of planes spanned by pairs of orthonormal basis elements |18j , the negativity of the Ricci 
scalar is only a sufficient (not necessary) condition for local instability of geodesic flow. For this reason, the negativity 
of the scalar provides a strong criterion of local instability. 



B. Information-Geometrodynamical Entropy 

At this point, we study the trajectories of the system on A4 S . We emphasize ED can be derived from a standard 
principle of least action (Maupertuis- Euler-Lagrange- Jacobi- type) p], [2l[ . The main differences are that the dynamics 
being considered here, namely Entropic Dynamics, is defined on a space of probability distributions M s , not on an 
ordinary linear space V and the standard coordinates of the system are replaced by statistical macrovariablcs M . 
The geodesic equations for the macrovariables of the Gaussian ED model are given by [22j , 

cPQv d&" dB p 

— -|- r^ 1 = (14) 

with /J, = 1, 2,...,6N. Observe that the geodesic equations are nonlinear, second order coupled ordinary differential 
equations. These equations describe a dynamics that is reversible and their solution is the trajectory between an 
initial and a final macrostate. The trajectory can be equally well traversed in both directions. We seek the explicit 
form of (Ti"4"l) for the pairs of statistical coordinates (jia , cri )■ Substituting the explicit expression of the Christof- 
fel connection coefficients into (fl4l) . the geodesic equations for the macrovariables and associated to the 
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microstate become, 

d?^ ] 2 d^ da^ _ <P*P 1 (da^V 1 ..,.„ 
dr 2 4*) dr rfr ' rfr 2 CT W I rfr J 2<r a) I dr > 

with a = 1, 2,...., iV and a = 1, 2, 3. This is a set of coupled ordinary differential equations, whose solutions are 




^ (T) = — ; — 2 \ (B <a) \ 2 ' (T) = — ; \ + c - (16) 

ex P (-2#>r) + ^ exp (-2#>r) + ^ 

The quantities -B Q a \ C« , (3^ are real integration constants and they can be evaluated once the boundary conditions 
are specified. We observe that since every geodesic is well-defined for all temporal parameters t, M. s constitutes a 
geodesically complete manifold [23}. ft is therefore a natural setting within which one may consider global questions 

and search for a weak criterion of chaos fl3| . Furthermore, since jia (j) < +oo and (r) < +oo Vt € M + , 

Va = 1, 2, 3 and Va = 1,.., N, the parameter space {©} (homeomorphic to A4 S ) is compact. The compactness of the 

configuration space manifold M B assures the folding mechanism of information-dynamical trajectories (the folding 
mechanism is a key-feature of true chaos, [HI). 

We are interested in investigating the stability of the trajectories of the ED model considered on Ai s . It is known 
[2~fl | that the Riemannian curvature of a manifold is closely connected with the behavior of the geodesies on it. If the 
Ricmannian curvature of a manifold is negative, geodesies (initially parallel) rapidly diverge from one another. For 

the sake of simplicity, we assume very special initial conditions: = A, (3^ = A e IR + , = 0, Va = 1, 2,...., N 
and a = 1, 2, 3. However, the conclusion we reach can be generalized to more arbitrary initial conditions. Recall that 

M. s is the space of probability distributions P (^X labeled by 6iV statistical parameters 0. These parameters arc 

the coordinates for the point P, and in these coordinates a volume element dVM s reads, 

N 3 p: 

dV Ms = ^gd eN = II II 7^T72^i Q) ^i tt) (17) 
a=i«=i 

where g = |dct The volume of an extended region AVm s (t; A) of Ai s is defined by, 

N 3 M< o) (r)o-(°)(r) 
tl.-l TT TT j I V- 



A^(r;A)^niI / / T^TI^ (18) 

where fia^ (r) and eri"^ (r) are given in (|16p and where the scalar A is the chosen quantity to define the one-parameter 
family of geodesies 3~g Ms (A) = f {©jtt s ( r ! ^)}a£r+' 67V - The quantity that encodes relevant information about the 
stability of neighboring volume elements is the the average volume Vm, ( t ! A) defined as (20| . 

T 

V Ms (r; A) = (AV Ms (r'; A)) T d = f 1 1 AVm. (r' ; A) dr' T ~°° e 3N ^. (19) 

o 

This asymptotic regime of diffusive evolution in (fT§|) describes the exponential increase of average volume elements 
on A4 S . The exponential instability characteristic of chaos forces the system to rapidly explore large areas (volumes) 
of the statistical manifolds. From equation (|19[) . we notice that the parameter A characterizes the exponential growth 
rate of average statistical volumes Vm s i T \ -M m -M. s . This suggests that A may play the same role ordinarily played by 
Lyapunov exponents. Indeed, it is interesting to note that this asymptotic behavior appears also in the conventional 
description of quantum chaos where the von Neumann entropy increases linearly at a rate determined by the Lyapunov 
exponents. The linear entropy increase as a quantum chaos criterion was introduced by Zurek and Paz [l5} . In our 
information-geometric approach a relevant variable that can be useful to study the degree of instability characterizing 
the ED model is the IGE quantity defined as [20], 

S Ms = Hm logV^ s (r; A). (20) 
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The IGE is intended to capture the temporal complexity (chaoticity) of ED theoretical models on curved statistical 
manifolds M s by considering the asymptotic temporal behaviors of the average statistical volumes occupied by the 
evolving macrovariables labelling points on A4 S . 
Substituting (fT5|) in fTO]) , equation (f2T))) becomes, 



S 



lim log 



/4, q) ( t 'K q) ( t ') 



nn 



V2 



'(0) 



'(0) 



(*») 



dr' 



WXt. 



(21) 



Before discussing the meaning of (|2ip . recall that in a rigorous examination of the entropy approach to the classical- 
quantum correspondence problem, Zurek and Paz consider the completely tractable model of an inverted harmonic 
oscillator with a potential V (x) = — n £ coupled to a high temperature (harmonic) bath. In their case, Q is analogous 
to a Lyapunov exponent in a genuinely chaotic system. On calculating the rate of change of von Neumann entropy, 
they show that [la ] 



. . t — >oo 
^von Neumann \J~ ) ~ 



(22) 



where S von Neumann (t) = —tr (p r (r) log p r (r)) is the von Neumann entropy of the system [24| and p r (r) is the reduced 
density matrix of the system at time r. The quantum entropy production rate is determined by the classical instability 
parameter Q. Given that the classical Lyapunov exponent to which f2 is analogous is equal to the Kolmogorov-Sinai 
(KS) entropy of the system [24| . this is indeed a remarkable characterization. It suggests that after a time, a quantum, 
classically chaotic system loses information to the environment at a rate determined entirely by the rate at which the 
classical system loses information as a result of its dynamics, namely, the KS entropy. 

In analogy to the Zurek-Paz quantum chaos criterion in its classical reversible limit (25|, we suggest a classical 
information-geometric criterion of linear IGE growth. The entropy- like quantity Sm 3 m (|21[) is the asymptotic limit 
of the natural logarithm of the statistical weight (AVm 3 ) t defined on A4 S and it grows linearly in time, a quantum 
feature of chaos. Indeed, equation (|21|1 may be considered the information-geometric analog of the Zurek-Paz chaos 
criterion. In our chaotic ED Gaussian model, the IGE production rate is determined by the information-geometric 
parameter A characterizing the exponential growth rate of average statistical volumes Vm s ( t ; A) in A4 S . 



C. Jacobi Field Intensity 

Finally, we consider the behavior of the one-parameter family of neighboring geodesies ^Fgm, W ^ 
{®M S ( r ; A )}AeR+" 6JV where > 

& (r; A) = - f - A3 , ,W (r; A) = (23) 
exp (-2At) + fa exp (-2Ar) + ^ 

with a = 1, 2,...., N and a = 1, 2, 3. The relative geodesic spread on a (non-maximally symmetric) curved manifold 
as A4 S is characterized by the Jacobi-Levi-Civita equation, the natural tool to tackle dynamical chaos (l8l. [26l|. 

K P ,^-se"^- = o (24) 



Dt 2 vp>? Q T Q T 

where the covariant derivative D D S ^T in (EH) is defined as [27| . 

Dt* ~ dr 2 + ^ dr dr + dr 2 + Q ^ dr dr ° + 

dB a dO* 3 

(25) 

and the Jacobi vector field J M is given by [22j , 

(t- \\\ 

SX. (26) 



def 



d®v (t; A) 
dX 
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Equation (|24|) forms a system of 6N coupled ordinary differential equations linear in the components of the deviation 
vector field (|26[) but nonlinear in derivatives of the metric ([7]). It describes the linearized geodesic flow: the lineariza- 
tion ignores the relative velocity of the geodesies. When the geodesies are neighboring but their relative velocity is 
arbitrary, the corresponding geodesic deviation equation is the so-called generalized Jacobi equation (28j . The non- 
linearity is due to the existence of velocity-dependent terms in the system. Neighboring geodesies accelerate relative 
to each other with a rate directly measured by the curvature tensor R a /3 7 s- Substituting (123|) in |24j) and neglecting 
the exponentially decaying terms in 50^ and its derivatives, integration of (|24p leads to the following asymptotic 
expression of the Jacobi vector field intensity, 

Jm s = \\J\\ = {9^J V Y T «°° %Ne XT . (27) 

Further details on the derivation of this result are in [2Qj. We conclude that the geodesic spread on M s is described 
by means of an exponentially divergent Jacobi vector field intensity Jm 3 , a classical feature of chaos. In our approach 
the quantity Aj, 

A j = f lim 

r — >oc 

would play the role of the conventional Lyapunov exponents. 

In conclusion, the main results of this work are encoded in the following equations, 

T — ^OO T — * CO \ 

TZ Ms = -3iV, S Ms ps 3NXt, J Ms ps 3Ne Xr . (29) 

The IGE grows linearly as a function of the number of Gaussian-distributed microstates of the system. This supports 
the fact that Sm s may be a useful measure of temporal complexity [2!|. Furthermore, these three indicators of chaotic- 
ity, the Ricci scalar curvature TZm^ the information-geometric entropy and the Jacobi vector field intensity Jms 

are proportional to 3A, the dimension of the microspace with microstates |a| underlying our chaotic ED Gaussian 

model. This proportionality leads to the conclusion that there is a substantial link among these information-geometric 
measures of chaoticity since they are all extensive functions of the dimensionality of the microspace underlying the 
macroscopic chaotic entropic dynamics, namely 

Ti-M s ~ $M S ~ J M s - (30) 

Equation (f3"0"]) represents the fundamental result of this work: curvature, information-geometric entropy and Jacobi 
field intensity are linked within our formalism. We are aware that our findings are reliable in the restrictive assumption 
of Gaussianity. However, we believe that with some additional technical machinery, more general conclusions can be 
achieved and this connection among indicators of chaoticity may be strengthened. 



log 



J Ms ( t ) 
J Ms (0) 



(28) 



IV. CONCLUSIONS 



In this paper, a Gaussian ED statistical model has been constructed on a 6iV-dimensional statistical manifold M. s . 
The macrocoordinates on the manifold are represented by the expectation values of microvariables associated with 
Gaussian distributions. The geometric structure of A4 S was studied in detail. It was shown that Ai s is a curved 
manifold of constant negative Ricci curvature —3N . The geodesies of the ED model are hyperbolic curves on Ai s . A 
study of the stability of geodesies on A4 S was presented. The notion of statistical volume elements was introduced to 

investigate the asymptotic behavior of a one-parameter family of neighboring volumes J 7 v Ms (A) = {Vm 3 (t; A)} AeR+ . 
An information-geometric analog of the Zurek-Paz chaos criterion was suggested. It was shown that the behavior of 
geodesies is characterized by exponential instability that leads to chaotic scenarios on the curved statistical manifold. 
These conclusions are supported by a study based on the geodesic deviation equations and on the asymptotic behavior 
of the Jacobi vector field intensity Jm 3 on M s . A Lyapunov exponent analog similar to that appearing in the 
Riemannian geometric approach to chaos was suggested as an indicator of chaoticity. On the basis of our analysis 
a relationship among an entropy-like quantity, chaoticity and curvature is proposed, suggesting to interpret the 
statistical curvature as a measure of the entropic dynamical chaoticity. We think this is a relevant result since a 
rigorous relation among curvature, Lyapunov exponents and Kolmogorov-Sinai entropy is still under investigation 
and since there does not exist a well defined unifying characterization of chaos in classical and quantum physics due 
to fundamental differences between the two theories. Finally we remark that based on the results obtained from 
the chosen ED model, it is not unreasonable to think that should the correct variables describing the true degrees 
of freedom of a physical system be identified, perhaps deeper insights into the foundations of models of physics and 
reasoning (and their relationship to each other) may be uncovered. 
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